o 
o 



o 

RANDOM MATRIX THEORY OF SCATTERING 
IN CHAOTIC AND DISORDERED MEDIA 

r/a ' Jean-Louis Pichard 



<D 



cd 



Service de Physique de I'Etat Condense, 
CEA-Saclay, 91191 Gif sur Yvette cedex, France 
pichard@drecam.saclay.cea.fr 



Published in: Waves and Imaging through Complex Media ed. by 
q . P. Sebbah, Kluwer Academic Publishers (2001). 

l cj; 

Abstract We review the random matrix theory describing elastic scattering through 
^. , zero-dimensional ballistic cavities (having chaotic classical dynamics) 

r-~^ ■ and quasi-one dimensional disordered systems. In zero dimension, gen- 

£C) ' eral symmetry considerations (flux conservation and time reversal sym- 

metry) are only considered, while the combination law of scatterers put 
C ' in series is taken into account in quasi-one dimension. Originally de- 

veloped for calculating the distribution of the electrical conductance of 
mesoscopic systems, this theory naturally reveals the universal behav- 
iors characterizing elastic scattering of various scalar waves. 

a 
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This chapter is a short introductory review of the random matrix 
descriptions of elastic scattering. Additional informations can be found 
q , in more exhaustive reviews [1, 2, 3]. The more recent reviews are given 

by Bohigas [4], by Beenakker [5] and by Guhr, Miiller-Groeling and 
Weidenmiiller [6] . The basic references for random matrix theory are the 
rS | book of Mehta [7] (see also Porter [8]) and the series of papers published 

C<3 ■ by Dyson [9] in 1962. 



1. GAUSSIAN ENSEMBLES OF HERMITIAN 

MATRICES 

For a statistical description of a matrix ensemble, one has to define a 
measure of the space of the matrices having the required symmetries. If 



one is interested by the distribution of a restricted set of parameters suit- 
able for describing the matrices (e. g. the eigenvalues of an Hermitian 
matrix), one has to use the system of coordinates using those parame- 
ters. The J acobian of the transformation (e. g. from matrix elements 
coordinates towards eigenvalue-eigenvector coordinates) yields correla- 
tions between those parameters. Those correlations (level repulsions) 
are at the origin of universal behaviors first observed in complex nuclei, 
then in small metallic particles, quantum billiards, hydrogen atom in 
a magnetic field, mesoscopic quantum systems, electro-magnetic cavi- 
ties... The simplest illustration is given by the Gaussian ensembles of 
Hermitian matrices introduced in this section. Another illustration is 
given in the following section: the distribution of the radial parameters 
characterizing a scattering matrix S or a transfer matrix M. 

For doing statistics with real numbers, one defines the probability 
P(dx) to have a real number x inside an infinitesimal interval of length 
dx: P(dx) = p(x)p(dx) where p(x) is a density and fi(dx) = dx the 
measure of an infinitesimal interval of the real axis. Similarly, for do- 
ing statistics with matrices X, one defines the measure p(dX) of an 
infinitesimal volume element dX of the matrix space in which X is de- 
fined and one gives a density probability p{X). The measure p(dX) 
is given by the symmetries of X, while the density p(X) may contain 
physical assumptions (e.g. minimum information density given a few 
physical constraints). 

For instance, let us introduce the Gaussian Orthogonal Ensemble 
(GOE) of real symmetric matrices H and the probability distribution 
of their eigenvalues Ei . A real symmetric matrix H = H = H* of size 
N has (iV 2 + N)/2 independent entries. The infinitesimal volume ele- 
ment dH has a measure p(dH) given by the product of the infinitesimal 
variations dHij of the N(N + l)/2 independent entries: 

N 

t i(dH) = l[dH ij . (1) 

A possible definition of the GOE density probability p(H) is given by a 
maximum entropy criterion. Minimizing [10] the information entropy 

I(p(H)) = -Jp(H)lnp(H)p(dH) (2) 

with an imposed expectation value for the trace of H 2 gives 

/ trH 2 \ 
p(#)ocexpl — 2^2"]- ^ 
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Since 

exp I — ^j- ) fjt(dH) = Y[ exp I -^-| J dH ti ]J exp I -— |- J dH^, 

the N(N + l)/2 independent matrix elements Hij (i < j) are uncorre- 
cted variables with Gaussian distributions of variance a 2 and c 2 /2 for 
the diagonal and the off-diagonal entries respectively. This ratio between 
the variances is important since it makes the GOE ensemble invariant 
under change of basis: p(H) is only function of the N eigenvalues Ei 
of H through J2i=i ^f an d does not depend on the eigenvectors of H. 
To calculate P(E±, . . . ,En), one has to go from the parameterization 
of H in terms of its matrix elements H^ to the parameterization of H 
in terms of its eigenvalue-eigenvector coordinates. The Jacobian of this 
change of coordinates is at the basis of the level repulsion and of the 
spectral rigidity characteristic of usual random matrix theories. 

A real symmetric matrix is diagonalizable by an orthogonal transfor- 
mation. Let us define the measure of an orthogonal transformation On- 
We first consider a 2d rotation O2 of angle 9. One has 

/ cos# sin# \ , . 

2 ~ \ -sin# cos6 ) W 

and by differentiation 

*» = ° 2 ( -m ) < 5 » 

Clearly, fi(d02) = d9 is the appropriate measure for the transforma- 
tion 0-2. The generalization to an arbitrary N x N orthogonal transfor- 
mation On is straightforward: 

N N = In (On + dO N ) T (0 N + dO N ) = In 
dO N = N dA N dA N = -dA N 
where In denotes the unit N x N matrix. dAN is a N x N real anti- 
symmetric matrix and /x(dOjv) = Ui<j dAij. 

H is diagonalizable by an orthogonal transformation O: H = OHdO t 
where Hjj is a real diagonal matrix of entries Ei and of measure [i(dH]j) = 
rii=i dEi. By differentiation, one gets 

dH = OdUO T (6) 

dU = dAH D -H D dA + dH D (7) 

where we have used dO = OdA and dA = —dA. The real symmet- 
ric matrix dH is related to dTC by an orthogonal transformation. The 



Jacobian is equal to one and /j,(dH) = T\i<jdHij = Yli<jd7iij. The 
product of the infinitesimal diagonal elements of dH. gives fi(dHi)), the 
off-diagonal contribution gives IIi<j \E% — Ej\dAy, and one eventually 
obtains the measure fi(dH) in terms of the measures n(dO) and ^(dHo) 

N 

n(dH) = Y[\Ei- E j \ f i(dH D ) f i(dO). (8) 

i<j 

In terms of the eigenvalue-eigenvector coordinates of H, the GOE dis- 
tribution becomes: 

P(dH)n(dH) = P(E 1 ,..., E N )n{dH D )p{dO) (9) 

where the joint probability distribution of the eigenvalues is identical to 
the Gibbs factor of a set of N point charges free to move on the real 
axis of the complex plane with a pairwise logarithmic repulsion and a 
quadratic confining potential at an inverse temperature (3 = 1: 



P(Ei,...,E N ) oc exp 



N N F 2 



(10) 



The pairwise repulsion coming from the Jacobian makes unlikely level 
degeneracies and explains why dramatically non-random an energy-level 
series really is. This is the Coulomb gas analogy usual in Random Ma- 
trix Theory. To appreciate how this random matrix approach is adapted 
to include symmetry breaking effects, let us assume that the matrix H 
is the Hamiltonian of an electron moving in a chaotic cavity. Applying 
a magnetic field removes time reversal symmetry, H becomes hermitian 
(H = H^) and fi(dH) = YlfLi dHiiTli<j dH^dHfj , taking into account 
the infinitesimal variations of the real and imaginary parts of its matrix 
elements (Hij = Hfj + iHfj). H is now diagonalizable by a unitary trans- 
formation U and dU = Uda where da is an infinitesimal anti-hermitian 
matrix (da = —da'), and /j,(dU) = Ili=i d a \ Ili<j dajjdafj. One obtains 
for hermitian matrices dH = UdTiU^ where dTi = daHn + dHr, — Hoda. 
The Jacobian of a unitary transformation being equal to one, one even- 
tually finds n(dH) = rij<j \Ei — Ej\ 2 /j,(dHr>) jx{dU) / (Oi dafi), the square 
coming from the fact that the non diagonal contribution of dfi is now 
complex. Breaking time reversal symmetry, one keeps the Coulomb gas 
analogy, with a temperature divided by a factor two (/5 = 1 — >■ 2). For 
electrons of spin 1/2, one can also break spin rotation symmetry (SRS) 
by spin orbit scattering, an effect which preserves time reversal sym- 
metry (TRS). The matrix elements of H are no longer real (/3 = 1) or 
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complex (/3 = 2), but quaternion real {(3 = 4) and the level distribution 
is still given by the Coulomb gas analogy with a temperature divided by 
a factor 4: (3 = 1 — ► 4. 

The main feature of those three Gaussian ensembles of random ma- 
trices is that p(H) does not couple eigenvalue and eigenvector variables. 
These ensembles are invariant under canonical transformations: orthog- 
onal transformations (f3 = 1) when the system invariant under TRS 
and SRS symmetries, unitary transformations {(3 = 2) in the absence 
of TRS and symplectic transformations {(3 = 4) with TRS and without 
SRS. The eigenvectors are totally random, the measure of the matrices 
O or U are given by the Haar measures over the orthogonal or unitary 
groups respectively, and the integration over the eigenvectors is trivial. 
This is the totally random character of the eigenvectors which makes 
the energy levels correlated (pairwise logarithmic repulsion) and subject 
to universal symmetry breaking effects (e.g. (3 = 1,4 — ► 2 when TRS is 
broken) . 

2. RADIAL PARAMETERIZATION OF 

SCATTERING MATRICES S AND 
MEASURES 

Similar random matrix theories can be adapted to matrices having 
different symmetries and can give the joint probability distribution of a 
subset of variables which can be used for their parameterization. An- 
other example is provided by the unitary matrices describing complex 
elastic scatterers. Let us consider a perfect waveguide characterized by 
iV quantized modes propagating to the right and N time reversed modes 
propagating to the left. Let us introduce in the middle of this wave guide 
a complex elastic scatterer described by its 2N x 2N scattering matrix 
S. This matrix describes the various transmission and reflection ampli- 
tudes to the right or to the left, has to be unitary for conserving the 
flux amplitudes (SS^ = S^S = I2N), and must be symmetric (S = S T ) 
if one has time reversal invariant scattering: 

S =( T t'r>) <»> 

The N x N matrices t and t' describe transmission amplitudes of the 
incoming fluxes to the right and left directions respectively, while r and 
r' describe reflections. Let us consider time reversal invariant scattering 
where S is symmetric and can be decomposed as 

S = U T .U (12) 

S + dS = U T (I 2N + idM)U (13) 




where U and dM are respectively a unitary and an infinitesimal real 
symmetric matrices. This decomposition is not unique, one can multiply 
U by an arbitrary orthogonal transformation O (U — > UO) but the 
measure /j,(dS) = IL<j dMy is uniquely defined since the Jacobian of 
the transformation dM — ► dM = OdMO T is equal to one. fi(dS) was 
expressed by Dyson [9] in terms of the eigenvectors and eigenvalues of S 
and their associated measures. The original motivation in Dyson's work 
was not at all to study a scattering problem, but to use the eigenvalue 
distribution of S for describing energy-level statistics (the eigenvalues of 
S being confined on the unit circle of the complex plane, one does not 
need the somewhat artificial GOE quadratic confining potential). For 
a scattering problem, the eigenvalue-eigenvector parameterization of S 
is not adapted and we introduce a more convenient one using 2 unitary 
N x N matrices u\ and u 2 and a diagonal N x N matrix A with iV 
real positive diagonal entries Ai, . . . , Ajv- Denoting T = (1 + A)" 1 and 
R = A(l + A)^ 1 , S can be written as: 

(14) 

In this parameterization, the transmission and reflection matrices be- 
come 

t = u 2 VTuJ t' = mVTuJ (15) 

r = —uivRu l r' = U2\Ru2 (16) 

and T and R contain the eigenvalues of tV and rr' (transmission and 
reflection eigenvalues). One can note that the transfer matrix M which 
gives the flux amplitudes of the right side of the scatterer in terms of 
the flux amplitudes of its left side can be written using the same param- 
eterization as 

M -{o u* 2 ){ Va VTFTa){ o u\ ) (17) 

M is pseudo unitary (flux conservation) 

and when one has TRS {(3 = 1), M must also satisfy the requirement: 
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M has the advantage to be multiplicative when one puts the scatterers 
in series. For the derivation of this parameterization of S, see Ref.[ll], 
and use the relation between the matrices M and S. To understand 
the interest of this parameterization, we introduce the variables Xi from 
Xi = sinh 2 X{ . One can see than M is decomposed in the product of a 
unitary transformation, followed by N hyperbolic rotations of angle xf. 

cosh Xi sinh X{ 
sinh X{ cosh Xi 

before a second unitary transformation. The N parameters \ are called 
the radial parameters of M (or S). 

Diagonalizing the matrix containing those parameters in the new pa- 
rameterization of S by the orthogonal matrix O 

n _ 1 ( {In-^/R) 1 ' 2 (In + ^R) 112 \ r2m 




and denoting 



and 



( ( ^ ^ ) (2.1.) 



U N 

-U N 



(22) 



one can write S = YY T where the unitary matrix Y = UOI. Defining 
the infinitesimal anti-hermitian and real antisymmetric matrices dA and 
dB from dU = UdA and dO = OdB, using dS = iYdMY T , one obtains: 

idM = dC + dC T (23) 

dC = I*(dB + T dAO)I (24) 

which allows us to write the Jacobian matrix of the change of coordinates 
dMij — ► (dAij , dBij ) . This matrix has a simple block diagonal form 
and its determinant gives the measure fj,(dS) in terms of the measures 
//(dA) = rij=i d\i and fi(dU) = Y\i = i^{dui). We have sketched the 
derivation when S is unitary symmetric {(3 = 1), but the extension 
to the three possible symmetry classes is straightforward and gives for 
Hp(dS) the general form [12] 

lip(dS) = PpiXu . . . , \ N )n{dk)n{dU) (25) 

P^Ai, . . . , \ N ) = exp -pHpiXu . . . , Ajv) (26) 



N N 

H p (X 1 ,...,X N ) = -Y / hi\X i -X j \ + ^2V l 3(X i ) (27) 

F /3 (A) = (AT + ^)ln(l + A) (28) 

If the scattering is not time reversal symmetric (/? = 2), 5 is no longer 
symmetric and one needs two additional unitary matrices U3 and 1x4 for 
parameterizing S and fi(dU) = IIj=i fJ>{dui). If one considers scattering 
of spin 1/2 particles by a TRS scatterer which removes SRS (spin-orbit 
scattering), (3 = 4. 

One can similarly show that the measure for the transfer matrices [2] 
is given in terms of the radial parameters A, by: 

N 2,4 

Hp{dM) = 1] \Xi - Xjfn(dA.) [] fi(dui) (29) 

i<j i=l 

3. ZERO DIMENSIONAL CHAOTIC 

SCATTERING 

Let us assume that the scatterer represented by the matrix S is a 
ballistic cavity of irregular shape having chaotic classical dynamics, as 
sketched in Fig. 1. We have in mind long scattering trajectories corre- 




Figure 1. Scheme of a cavity giving rise to zero-dimensional chaotic scaterring. 

sponding to particles reflected many times inside the cavity before being 
transmitted or reflected. Let us assume that the shape of the cavity is 



Random Matrix Theory of Scatteringin Chaotic and Disordered Media 9 

slightly modified, or that the wave length of the incoming fluxes varies. 
The scattering will be deeply re-organized and a statistical description of 
the fluctuations of the scattering amplitudes becomes necessary. To this 
end, we need to define a statistical ensemble of scattering matrices S and 
we will assume that the scatterer will visit this ensemble when one varies 
a tunable parameter (shape of the cavity, wave-length, applied magnetic 
field reorganizing the quantum interferences if we consider electron elas- 
tic scattering). The simplest ensemble is the one where all the scatter- 
ing processes are equiprobable, which does not contain any information 
about the system excepted its basic symmetries. Those ensembles of 
minimum information entropy for S are the circular ensembles [9] in- 
troduced by Dyson in 1962, for which the probability to find S inside a 
volume element dS is 

P(dS) = ^Hp{dS) (30) 

where Vp is a normalization constant. One obtains for the radial pa- 
rameters Aj a Coulomb gas analogy very similar to the GOE-GUE-GSE 
Coulomb gas analogies for the energy level of a random Hamiltonian, ex- 
cepted two noticeable differences: (i) the \ are real positive in contrast 
to the Ei which are only real, the Coulomb gas is free to move only on 
the positive part of the real axis in the complex plane (ii) the confining 
potential Vg(A) is implied by the symmetries of S (in contrast to the 
quadratic potential given by a certain choice of p(H)) and depends on 
the symmetry parameter /3. 

Let us see the implication for the total transmission probability T = 
trtv = J2i=i Ti = Z)i=i(l + ^i) -1 - When S gives Fermi wave scattering 
by a mesoscopic scatterer coupling to electron reservoirs, T gives its 
electrical conductance in units of 2e 2 /h (Landauer formula). 

If N = 1 (single mode wave guide), the probability distribution of T 
exhibits strong symmetry breaking effects (see Fig. 2): 

P(T) = ^T 3 / 2 - 1 (31) 

If one measures the electrical conductance of a chaotic cavity (quan- 
tum billiard) coupled to two electron reservoirs by single mode leads, 
the conductance is likely to be close to if one has TRS and SRS, to 
have a value uniformly distributed between and 2e 2 /h if one applies a 
small magnetic field (no TRS), and a value close to 2e 2 /h if the electron 
reflection on the walls of the cavity is accompanied by spin-orbit scat- 
tering (no SRS). The distribution P{T) can be explored by changing the 
Fermi energy with a metallic gate, by small deformations of the cavity 
or by applying a magnetic field. 
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Figure 2. Symmetry breaking effect on the distribution P(T) for a chaotic cavity 
coupied to leads via single mode contacts (N = 1). 

If N >> 1 the symmetry breaking effects are much smaller. A first 
one is the suppression, when there is no TRS, of a small "weak localiza- 
tion" correction to the average transmission < T >, the second effect of 
removing TRS is to halve the universal variance of T. 

3.1. WEAK LOCALIZATION CORRECTIONS 

The large N limit of the density p{\) 

N . . N 

PpW = T, •••/ I[d\ i Pp(\ 1 ,...,\ N )6{\-\ i ) (32) 

can be calculated using an equation derived by Dyson [13] 

pp{\')d\i (i-2d\npp{\) = dV p {\) 
R + A -A' 2/3 d\ dX l ' 

Taking the V^(A) characterizing the three circular ensembles for S, one 
gets for the ensemble averaged total transmission: 

The term oc A^ is obvious: having chaotic scattering, the probability to 
be reflected equals the probability to be transmitted. However there is 
a small correction of order 1 which reduces by —1/4 T when there is 
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TRS and SRS, which disappears without TRS and which enhances T 
by a factor 1/8 without SRS. This is the analog in a quantum billiard 
of the well-known weak-localization corrections to the Boltzmann-Drude 
conductance of a disordered system. 

3.2. UNIVERSAL CONDUCTANCE 

FLUCTUATIONS 

One of the main phenomena which is naturally explained by random 
matrix theory [14] are the "universal conductance fluctuations" (UCF) 
characterizing mesoscopic conductors. When one varies the "conduc- 
tance" T with an external parameter (Fermi wave length, magnetic field 
...) one generates fluctuations of magnitude independent of < T >. To 
calculate < 5T 2 > one needs to know the density-density correlation 
function of the A-parameters. In the limit N — > oo, one can simply cal- 
culate [15] this variance. Exploiting the Coulomb gas analogy, one can 
write 

K p (X, A') =< J2 5 ( X ~ X i) S ( X ~ X j) > -PpWPpM (35) 

ij 

as a functional derivative: 

When N — ► oo, Vp(X) — ► J R+ ln|A — X'\pp(X) + const (this amounts to 
neglect the term responsible for the weak localization correction in Eq. 
(33). In the large TV-limit, Vp{\) becomes a linear functional of p(X'). 
This implies an important result: Kp(X, A') depends only on the nature 
of the pairwise repulsion between the radial parameters and becomes 
independent of the confining potential Vp(X). The evaluation of the 
functional derivative giving Kn(X,X') is straightforward: 

lim Kp(\,\') = Lln ^~^ . (37) 

iV-oo ^ 7T 2 /3 VA + VA 

to eventually give 

*2,„w . - 1 f f ;x .„/,.. V^-vtfda(A)da(A') 



<5\A)>^— 1 / dXdX'ln " v \> \/ (38) 

(3tt z Jr+ Jr+ ^/x + ^/\ dX dX' 

for the variance < 5 2 A > of a linear statistics A = J2i=i a (^i) of the 
A-parameter. Taking a(X) = (1 + A) -1 , one gets 
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The variance of 7" is just a number which depends on [3 but not on 
<T>. 

When N is finite, weak localization corrections and variances can 
be exactly calculated using a method introduced by Mehta, Gaudin 
and Dyson. Using a set of orthogonal polynomials {(3 = 2) or shew- 
orthogonal polynomials {[3 = 1,4), one can perform the integration of 
P^(Ai, . . . , Ajv) over an arbitrary set of A-parameters in order to have 
Pb(X), Kg(X,X') and higher order correlation functions. Going to the 
variables Tj = (1 + Aj) the polynomials required [12] to perform the 
integrals for [3 = 2 are the Legendre polynomials. 

4. MANY CHANNEL DISORDERED WAVE 

GUIDE 

So far, we have done the simplest random matrix exercise where all 
the scatterers are taken with a uniform distribution. Those circular 
ensembles are suitable to describe scattering in zero-dimensional chaotic 
cavities, as it has been numerically checked. Another exactly solvable 
case is provided by a Id series of scatterers. Let us consider a quasi- Id 
disordered wire (or wave guide) with N channels, as sketched in Fig. 3. 
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Figure 3. Scheme of a quasi- Id disoredered wave guide of length L » L t with 
N = (Lt/Xpf- 1 modes. 

Typically, if the (Fermi) wavelength is A^, and L t ~ the transverse 
section, one has N = (L$/Air) d_1 . Let us consider a "building block" of 
length 6L and of transfer matrix M^l ■ A suitable statistical ensemble for 
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M§l is given by a maximum entropy ensemble where < T > is imposed. 
If 5L is larger than the elastic mean free path I, a natural requirement 
is to impose Ohm's law (the conductance < g >oc< T >oc Nl^dL)^ 1 ). 
This gives for the building block P(M$ L ) oc exp-(A < T >)/j,(dM), 
where A is the Lagrange multiplier associated to the imposed constraint. 
This ensemble preserves the logarithmic pairwise interaction between 
the A-parameter and the distribution of the auxiliary unitary matri- 
ces Ui derived for the circular ensembles, but gives rise to a differ- 
ent potential Vg(X). Considering a Id series of such building blocks, 
of length L/l, and exploiting the multiplicative combination law of M 
(Ml+&l = M^.Msl), one can derive a Fokker-Planck equation [5, 16, 17] 
for Pp(\i,...,\ N ,L/l): 



fj V -J. 7 7 

dL 


■^ = DA x P p (\ l ,...,\ N ) 
D 2 


(40) 
(41) 

(42) 
(43) 


A A : 


" {(3N + 2 - (3)1 

= E^(i + w^- 7 

N 

j = H\X i -\ j f 



i<j 

In this statistical description of M (or S), the radial parameters and the 
matrices Ui are statistically uncorrelated. The U{ remain distributed with 
Haar measure over the N x N unitary group. Increasing the number L/l 
of blocks put in series only changes the statistics of the radial parameters 
as given by the above Fokker-Planck equation, This is the limitation 
of this "isotropic" model which allows us to solve it entirely and to 
describe quasi-lci localization. But the transverse system dimensions 
appear only through the parameter TV". A strip of purely Id transverse 
section and a bar with a 2d section are treated the same way: within 
the Od approximation for the transverse dynamics. When N = 1, J = 
1 and the resulting equation was originally derived [18] in 1959 in a 
work entitled: "waveguides with random inhomogeneities and Brownian 
motion in the Lobachevsky plane". The nature of the Fokker-Planck 
equation was indeed correctly identified as a "Heat equation in a space 
of negative curvature" . Looking in the mathematical literature, one can 
realize that A\ is the radial part of the Laplace-Beltrami operator in 
the space of the transfer matrices M. Increasing L yields for M an 
"isotropic" (invariant under the unitary transformation Ui) Brownian 
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motion in the transfer matrix space which is only characterized by the 
diffusion constant D. 

4.1. QUASI-1D LOCALIZATION 

When L increases and exceeds the localization length £, the system 
becomes an Anderson insulator and it is more convenient to use the 
variables Xj = L/& (\ = sinh x{), where the lengths £j characterize the 
exponential decays of the transmission channels of the quasi-ld scatterer 
(T-i rs exp — (2L/£j)). When L — ► oo the variables l/£j are given by the 
Lyapounov exponents of the multiplicative transfer matrix M, giving 
N decay lengths £jv << 6v-i << £i, the largest of them defining [19] 
the localization length £ of the quasi- Id system. The Fokker-Planck 
equation becomes: 

? ^r^-'>iU p ^O (44> 

N -. N 

Q = yj In | sinh x i — sinh Xj\ — — \] hi | sinh Xj | (45) 

i<j i=\ 

where the different channels are coupled via O. However, when L — > oo, 
| sinh xi — sinh Xj | ~ exp 2xj when Xj > > Xj , the channels become 
decoupled and fi -► -(2//3) E^l + PN - (3)xi. The Fokker-Planck 
equation becomes solvable and gives [1]: 



Pf3(xi,...,X N ,L) = ( 

where & = (70/(1 + /9JV - 
important results: 

(i) a universal symmetry breaking effect for the localization length £ 
when (3 = 1 — ► 2: 

£ = (/?iV + 2 - /?)/ (47) 

(ii) a normal distribution for In 7" where —2 < In 7" >=< <5 2 ln7~ >. 
Those symmetry breaking effects have been observed in magneto-transport 
measurements performed in disordered wires [20, 21, 22], where the con- 
ductance g oc exp -(2L/£) and £= (N + 1)1 -> 2NI when f3 = 1 -> 2. 
However this theory neglecting electron-electron interactions, it would 
be important to check those universal symmetry breaking effects in the 
localized regime using other waves (light, sound ...). 



ll ,n y-r ( l l , 
1=1 


- If) <«> 


- /?) and 7 = 0N + 2 - 


- /?. One gets two 
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4.2. MAPPING ONTO A 

CALOGERO-SUTHERLAND MODEL OF 
INTERACTING FERMIONS 

When L < £, the system is a disordered conductor where a number 
N e ff < N of channels are still opened, giving < T >~ Nl/L (Ohm's 
law) if one ignores the (small) weak-localization corrections. Approxi- 
mations can be used also in this limit (L « £) to calculate the quasi-l(i 
weak-localization corrections ST = {(3 — 2)/ (3/3) and the UCF variance 
2/(15/3). One can notice that those values valid for the quasi-ld wire 
are close, but not identical to those derived from the circular ensem- 
bles. The small difference between the UCF variances tells us (see the 
arguments given for having K(X, A') in the large iV-limit) that the radial 
parameters cannot have [15] exactly the pairwise logarithmic repulsion 
in quasi- Id. This can be understood when (3 = 2 where one can solve 
the Fokker-Planck equation for any values of N, using a transformation 
originally introduced by Sutherland to solve Dyson's Brownian motion 
model. The distribution Pp(xi,...,XN,L) is related [5, 23] to a wave 
function \P(a;i, . . . , xn, L) by the transformation 

P = *exp-(^) (48) 

and the Fokker-Planck equation for P becomes a Schrodinger equation 
for ^ in imaginary time. 

~ l dl = H ^ (49) 

where the Hamiltonian is given by 



H = — El— + ■ , 2 „ )+(0-2)U(x i ,x j ) (50) 

, . , sinh 2 2xj + sinh 2 2x, , rn . 

U(xi,Xj) = — > g ^7 (51) 

27 - 1 . (cosh2xj — cosh 2xi) 2 

For (3 = 2, the "particles" do not interact and the equation can be 
solved to give [23] a small change in the Coulomb gas analogy for the 
radial parameters: 

1 ,1 2 ,— 2 1— 

— In \\j — A, — ► — In A,- — A, In arsinh vA; — arsinh \ \\\ 

I l J \ cy \ <■ Jlol VI V 

2t ^_ 2 , 2 




yielding the change < S T >= y4g — v y#g. 
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For (3 7^ 2, one has not yet found how to solve the Schrodinger equa- 
tion. Let us note that the Calogero-Sutherland Hamiltonian derived 
from the original Dyson's Brownian motion model is now solved [24] 
when (3/2 = p/q where p and q are integer. 

5. SUMMARY 

We have seen how to derive the probability distributions of S and how 
to extract the distribution of T using very little information (symmetries, 
combination law in quasilti). This can be simply done as far as the radial 
parameters are decoupled from the auxiliary matrices Ui, which limits 
the method to Od and quasi-ld A more difficult task remaining to be 
achieved is to go beyond this limit, and to describe elastic scattering 
in 2d and 3d, including possible Anderson localization. One can ask to 
what extend a real chaotic cavity or a quasi-ld disordered wave-guide 
is accurately described by those random matrix theories. Numerical 
calculations [25] of the distribution of T in (suitably) designed chaotic 
cavities (suitably) connected to two many-channel ballistic waveguides 
confirm the random matrix results. For quasi- Id disordered wires, a field 
theory approach has been derived [26] assuming local diffusive dynamics. 
One obtains a non-linear c-model which gives using supersymmetry < 
T(C) > and < 5 2 T(L) > in the large iV-limit. The results turn out to 
coincide [27, 28] with those given by the Fokker-Planck equation in this 
limit. The finite N behaviors which are calculable by the Fokker-Planck 
equation are still out of reach using the cr-model, as well as the moments 
of InT which are the statistically meaningful observables (related to a 
normal distribution) in the localized limit. 
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